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Abstract 



The fully-packed loop model of closed paths covering the honeycomb lattice is 
studied through its identification with the sl q (3) integrable lattice model. Some known 
results from the Bethe ansatz solution of this model are reviewed. The free energy, 
correlation length, and the ensemble average loop length are given explicitly for the 
many-loop phase. The results are compared with the known result for the model's 

o 

O ■ surface tension. A perturbative formalism is introduced and used to verify results. 
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Figure 1: An example of a fully-packed loop configuration. With periodic boundary condi- 
tions, this configuration has 11 loops. 

1 Review of the FPL model 

In a recent article PJ] Blote and Nienhuis performed numerical investigations of what they 
termed the fully-packed loop (FPL) model. This is a statistical model where the ensemble 
is the set of all combinations of closed paths on the honeycomb lattice that visit every 
vertex and do not intersect. The Boltzmann weight of such a filling set of paths is just the 
exponential of the number of paths, i.e. the energy of a configuration is the number of closed 
loops used to cover the lattice. An example of a fully-packed configuration of loops on this 
lattice is shown in Figure [p. The partition function for this model may be represented as 



Z FPL {n) = Y^ 



n 



P(C) 



XD 
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where the sum is over all C, the coverings of the vertices of the hexagonal lattice by closed 
nonintersecting paths, P(C) is the number of paths in the covering C, and n is a generalized 
activity. This model was originally studied for its interest as the low-temperature limit of 
the 0(n) vector lattice models |2|, HJ. In this limit, the dimensionality of vectors n is just 
the activity n in equation ( |1 . 1|) . 

The partition function ( |1.1|) is apparently the generating function for the numbers of 
ways to cover the hexagonal lattice by any number of closed paths. Its calculation in the 
thermodynamic limit is an interesting combinatorial problem. 

More recently, Batchelor, Suzuki and Yung M pointed out that previous authors |5]. 
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Figure 2: The definition of vertex weights in the integrable lattice model associated to 
sl q (d + 1). 7 and 9 are parameters of the model. States on edges are labelled by roman 
indices with (d + 1) possible values. In the formula, A a & = r^fr + 2^rf . 

had exploited an identification of the FPL model with the integrable lattice model associated 
to the quantum group sl q (3) [0. This integrable model is a vertex model on the square lattice 
where each link of the lattice can be in one of three states, and the vertex weights are given 
by the R-matrix for sl q (3) as in Figure §. The R- matrix depends on a deformation parameter 
7 = log(g), as well as a spectral parameter 9 typical of integrable theories. Denoting the 
sl q (3) partition function as Z s i q ^(j,9), the precise idenfication is 



Z FPL (e^ + e-T) = {ei - e^y N Z slq (3) ( 7 , 9) 



1.2) 



where iV is the volume of the lattice (the number of hexagonal faces). Since the model 
is integrable, much exact information can be derived. In particular, the model's Bethe 
equations have been constructed and solved. 

One of the more important results that have been derived in this way is the existence of 
a phase transition in the model ( |1.1|) at n — 2 || |8|]. At larger n, larger numbers of loops are 



favored and at smaller n configurations with fewer loops are favored. It has been conjectured 
that this transition is between a large-n phase where the average loop length is finite and a 
small-n phase where this average is infinite. 

In section |] a simple relation between the free energy of the FPL model and the ensemble 
average length of loops is derived. From the known solution to the Bethe equations of the 
sl q (3) integrable lattice model, the free energy is identified and used to graph the exact value 
of the average loop length as a function of n. 

Identifying n = e 7 + e~ 7 , it becomes apparent that n > 2 corresponds to the integrable 



model for 7 real and < n < 2 corresponds to 7 purely imaginary. The former phase is 
known [0] to be massive, in the sense that there is a gap in the spectrum of eigenvalues 
of the transfer matrix between the leading eigenvalue and the next-leading eigenvalue. By 



standard arguments [|10[ , this implies a finite correlation length. The gap tends to zero as 7 
goes to zero, showing that n = 2 is a critical point of the model (|1.1| ). 

In section [| this correlation length is studied by considering the spectrum of eigenvalues 
of this transfer matrix. The spectrum may be deduced directly from the model's Bethe 
equations. In the case that the transfer matrix is symmetric and therefore has real eigen- 
values, the correlation length is related to the maximum eigenvalue A max and next-leading 
eigenvalue A x by 

r^iog^f. (1.3) 

2 Review of sl q integrable models 

The R-matrix of the sl q (d + 1) quantum group is an (n + l) 2 x (n + l) 2 matrix that may 

be interpreted as a matrix of Boltzmann weights of a vertex model on the square lattice 

as shown in Figure |2|. Since this matrix satisfies the Yang-Baxter equation, the associated 

transfer matrix commutes with itself evaluated at differing values of the spectral parameter 

and the model is exactly solvable by a recursive set of d nested Bethe ansatze 0. 

The formula for eigenvalues of the transfer matrix for the range of parameters, 9 > — 17 

and 7 real and positive, is 

a _ TT sinh ^ Afc + h~ d ) (9U 

A - J, 1 sinh(iA fe - §7 - 9) [ } 

where the product is over roots, A& of a set of Bethe equations and we have neglected terms 
that do not contribute in the thermodynamic limit. In this limit the A& are distributed in 
the interval [— 7r/2, it/ 2] with the density, 

p„„(A) = t e 2 - Ai • Sgrr i - < 2 - 2 > 

m f^oo TT smh[(d + 1)7717] 

For eigenvalues near the maximum eigenvalue, the changes in the distribution are parame- 
terized by the locations of holes 9\, q = 1 • • • d, h = 1 • • • N q , according to 



7r sinh[(d + 1)7777] 
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The numbers of holes are constrained to satisfy the relation, 

d 

J2 qN q = multiple of (n + 1). (2.4) 

The formulas ( |2.1| ), ( p.3| ) may be combined in the thermodynamic limit to yield a formula 
for eigenvalues of the transfer matrix of the integrable model, 



logA=/ p(A)lo. 
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3 Free energy and average loop length 

The formulas of the preceding section in the case d = 2 yield directly the free energy density 
of the model ( |1.1| ) for the n > 2 phase as the logarithm of the maximum eigenvalue of 
the transfer matrix, rescaled by the factor of equation ( |1.2|) . This free energy was actually 
derived in 1970 by Baxter |§ as the solution to a weighted three-coloring problem on the 
honeycomb lattice. The free energy density of the FPL model in the n > 2 phase is 

F F PL{n) = lim — log Z FPL (n) 

f °° n - «-6m+2\2 ^1 

= 4^ (i-UV-U } (3 - i) 

where n = q + q' 1 , and q = e 7 > 1. This function has an essential singularity at q — 1. The 
free energy for n < 2 and with periodic boundary conditions is given in integral form in [|J. 
It is interesting to note that for both phases, the free energy density gives the ensemble 
average length of loops. Since a configuration C on a lattice of iV faces has 2N occupied 
links, the total length of loops is always 2N . The average loop length of configuration C is 
therefore 2N/P(C). If we define the ensemble average loop length L^in) by 



then from inspection of equation ( |1 . 1|) it is clear that 

— [L N (n)Z FPL {n)} = Z FPL {n). (3.3) 

an n 

The general solution to this equation can be written up to quadrature by direct integration: 

1 r n 27V 

Ln(u) = / —Z FPL (n')dn' (3.4) 
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Figure 3: This is a graph of the average loop length of the FPL model versus n, the fugacity 
of loops. The critical point is at n = 2. 



where the lower limit of integration is an undetermined constant. In terms of the free energy 
density Ffpl = 0-/N) log Zfpl, this becomes 

,-n e NF FPL (n') 
IC 



r n p NF FPL {n') 

L N (n) = 2Ne- NFFPL{n) / dn'. 

Jc n 



(3.5) 



The integral in equation ( |3.5[ ) can be evaluated by steepest descent. The result is 

exp 



t n 2N 
J-'N\P') — exp 
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The constant of integration may now be determined from the known value of L^ip) at 
n — > oo. As will be shown in section |5], in this limit Z^(n) ~ 3n N ^ 3 , FPL (n) ~ l/3n, and 
Ln{it) — 6. These imply that C = — oo, so in the thermodynamic limit 

2 



L N (n) 
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dF F 



(3.7) 



dn 



In this calculation we have neglected corrections of order 1/7Y to L^{n). 

A graph of the ensemble average loop length versus n in the large-n phase is shown in 
Figure 0. This verifies the conjecture of B that loop length diverges at the critical point. 



4 Correlation length 



To obtain the correlation length, we must compute the expression (|2.5| ) for the minimal hole 
distribution. When d = 2, there are two choices for the N q . Either Ni = 3 and N 2 = 0, or 
Ni = 1 and A 2 = 1. In each case, the eigenvalue gap is minimized for holes at d\ = | where 
the sum in equation (|2.3| ) after integration in Q2.5|) is oscillatory. The transfer matrix of the 
model is symmetric at the point 9 = —\ r ) 1 and its eigenvalues are then real. After setting 9 
to this value, the correlation length of the model is given by equation ( |1.3j ). 

Considering the case iVi = 3 and N% = 0, we denote the next-leading eigenvalue for this 
hole distribution as A30. The equation (|2.5|) together with the formula for densities (pT 
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The integral in equation 



introducing the variables q = e 

sinh(7A + 57 



2|) can easily be performed by contour integration. After 
7 and z = e 8 , the result for — ^7 < 9 < is 
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Substituting this result into equation 
A30 



log 
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m>0 '" \Q 1 

After expanding the demoninator of the summand in a power series in q~ x , this may be 
resummed to the form, 
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This form is now convergent at the symmetric point, 9 = — ^7 or equivalently z 2 = q l . We 
may therefore evaluate it there to obtain the correlation length according to equation ( |1.3| ), 

'l + q 2 q- Qm ){l + q A q- &m )\ 
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This is the desired result, the correlation length of the FPL model where n = q + q l and 



q > 1, or equivalently q = +\ / n 2 — 4. 

The other possible choice of holes, N\ — 1 and N2 = 1 may be computed in the same 
way to give 

Amax _ . f (1 + g 2 g-6")(l + g 3 g- 6m )(l + q 3 q~ 6m )(l + q 4 q~ 6m ) ] ,. ?) 

° g An ° g rii (l + g- 6m )(l + gg- 6m )(l + gV 6m )(l + gV 6m ) J" l ' J 

This quantity is greater than (|4.6|) for all q > 1, so it is not the inverse correlation length. 



For large g, the inequality may be seen by considering the limiting forms of expressions (4.6) 
and ( |4.7|) . Rigorously, the multiplicands in ( |4.7| ) may be seen to be greater than those in 
(|4.6|) term by term in m. 



5 Perturbative analysis 

The FPL model has a natural large-n expansion which allows simple perturbative verifica- 
tions of results. 

When n is large, the dominant configurations are those with large numbers of loops. The 
perturbative procedure is to approximate the sum over states by including the configurations 
with the highest numbers of loops. 

On a hexagonal lattice with number of faces N a multiple of three, there are three 
configurations with the maximum possible number of loops. In these states, one out of every 
three faces has a small loop around it and these small loops lie on a triangular lattice. A 
sample is shown in Figure |j. These three configurations differ by translations and each has 
iV/3 loops. 

The smallest change in the number of loops that can be made is to introduce a defect 
somewhere in one of the maximal configurations, as shown in Figure [|. There are 2N/3 
different such defects that can be introduced and each reduces the number of loops by 2. 

Introducing defects in this way, we can reach all possible configurations. To see that 
this is so, we can represent a configuration by labelling the links on the lattice that do not 
contain part of a path. One of every three links is unoccupied, and every vertex touches one 
unoccupied link. These unoccupied links form a dimer configuration for the vertices of the 
lattice. If we draw rhombuses around every dimer and interpret the resulting picture as the 
projection of the edges of a stack of cubes, we see that a FPL configuration is equivalent to 
a stack of cubes. Such an identification is shown in Figure || 
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Figure 4: A sample from a configuration with the maximum number of loops. 
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Figure 5: A configuration with two fewer than the maximum number of loops. 




Figure 6: An example of the identification of FPL configurations and stacks of cubes. One 
rhombus is drawn centered on each unoccupied link. 



In this new representation, the action of inserting a defect is just the action of adding or 
removing a cube. This identification is exhibited in Figure [7[ The result then follows that 
since every stack of cubes can be made by adding or removing cubes, every FPL configuration 
can be made from one of the maximal ones by inserting some combination of defects. 

To obtain an approximation for the free energy, consider first the maximal state shown 
in Figure |]. For a lattice of iV faces, this configuration has (N/3) loops. There are 3 
such configurations corresponding to the three-fold translational degeneracy of the state. To 
lowest order then Zfpl = 3n N ' 3 [l + 0(n~ 1 )]. This result was used in section [3] to determine 
the asymptotics of the average loop length. 

Allowing defects, there are (2N/3) locations for a defect and each defect reduces the 
number of loops by two. Defects may be applied in any number and in any combination, 
giving the usual sum over disconnected diagrams. We can write this as the exponential of 
the connected diagram (one defect) and we will be correct except for the effects of excluded 
volumes which begin with two-defect connected diagrams and are therefore higher order. To 
the next order, Z FPL = 3n N / 3 exp[(2N/3)n- 2 ] exp[0(n~ 4 )]. 

Perturbatively calculating the FPL free energy, we see that 



FfplM 



log(n) 



2n 



+ 0(n 



3 ov ' 3 

in conformity with Baxter's result shown in equation ( [3.1|) . 
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Figure 7: In the cube representation, introducing a defect is adding or removing a cube. 



This point of view incidentally leads to a simple expression for the entropy density of the 
FPL configurations at n — 1. At this point, all configurations are weighted equally and Zfpl 
is just the number of configurations, or the exponential of the entropy. Then calculating the 
partition function is just the problem of counting the number of coverings of the honeycomb 
lattice by paths, which is the number of different possible stacks of cubes, which is the old 
combinatorial problem of counting plane partitions. Elser |12| has calculated the asymptotics 
of plane partitions for large arrays of numbers. 

The result applied to this case is entirely dependent on the shape of the boundary, even in 
the thermodynamic limit. This is to be expected when n = 1, because this is in the small-n 
phase where the model is critical. For a lattice of N faces and free boundary conditions, the 
maximum entropy is obtained for a hexagon-shaped boundary and in that case the partition 
function is asymptotically 

'3 



ZfpM = exp 



N -log3-21og2 



(5.2) 



11 
6 Comparison with surface tension 

The ground state of the sl q (d) integrable lattice model is (d+l)-fold degenerate. This implies 
the existence of a notion of interfacial tension £(7) away from the critical point between 
regions of differing antiferromagnetic polarization. By considering finite-size corrections, 
de Vega || has derived transcendental equations for this interfacial tension and computed 
the asymptotic behavior of S in the limits 7 — > and 7 — > 00. 

Scaling arguments originally due to Widom [|ll|] predict that the scaling relation, S£ ~ 1 
should hold near the critical point, 7 — > or equivalently n — > 2 + . It would be interesting 
to test this relation in this case, but we know of no explicit expression for the interfacial 
tension. 

Away from the critical point however, a comparison can be made. The asymptotic 
behavior of the interfacial tension for 7 — > 00 was extracted by deVega, and the result is 

S( 7 ) = ^ T 7 + 0(1). (6.1) 

In the case of the FPL model, d = 2, n = e 7 + e -7 , and 

S(n) = hog(n) + 0(l). (6.2) 

This result may be compared with a perturbative calculation. Consider the sum over FPL 
states at large-n with the constraint that boundary conditions are fixed to cause frustration 
in the bulk, as in Figure ||. The configuration in that figure has the maximum number of 
loops possible and is the analog of the configuration shown in Figure |j. Denoting the sum 
over defects in this configuration by Z' FPL , the interfacial tension is defined to be the change 
in free energy per unit length of the interface: 

f^ ~ e~ LS (6.3) 

%FPL 

where L is the vertical size of the lattice. 

For a lattice of N faces, the maximum number of loops possible in the presence of the 
constraint is (N/3) — (2L/3) instead of (N/3). The maximal state in the presence of the 
constraint is now 3 x 2 2L / 3 -fold degenerate, because there are (2L/3) locations near the 
interface where defects may be freely introduced without changing the number of loops. To 
lowest order therefore, 

Z FPL ~ 3n N/3 (6.4) 

Z' FpL ~ 32 ^/3 n (N-2L)/3_ (65) 
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Figure 8: An interface separating two regions of differing polarization. 



Reading off the exponents, we have from equation (|6.3|) the result that 

S( 7 ) = |bg(n)+0(1). (6.6) 

Equation ( |6.6| ) is apparently consistent with the large-7 asymptotics derived in ||. 

Equation ( |6.6|) together with equation ( |4.6|) show that S£, 7^ 1 in the FPL model. More 
generally, from the correlation length calculation it is clear that for large 7 the leading 
behavior of the correlation length for any value of d will always be 7, and the leading 
behavior of S is always 'yd/ (d + 1). 
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